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We study minimal surfaces which arise in wetting and capillarity phenomena. Using conformal coordinates, 
we reduce the problem to a set of coupled boundary equations for the contact line of the fluid surface, and then 
derive simple diagrammatic rules to calculate the non-linear corrections to the Joanny-de Gennes energy. We 
argue that perturbation theory is quasi-local, i.e. that all geometric length scales of the fluid container decouple 
from the short-wavelength deformations of the contact line. This is illustrated by a calculation of the linearized 
interaction between contact lines on two opposite parallel walls. We present a simple algorithm to compute the 
minimal surface and its energy based on these ideas. We also point out the intriguing singularities that arise in 
the Legendre transformation from the pure Dirichlet to the mixed Dirichlet-Neumann problem. 



I. INTRODUCTION 



Minimal surfaces, i.e. surfaces of minimal area with specified boundary conditions, are found in many areas of physics, 
. mathematics and biology. Their existence, uniqueness and other properties (such as possible singularities or stability) are still 
actively studied by mathematicians jT|]. In the laboratory, minimal surfaces are most commonly realized as soap films bounded 
by a given wire-frame, a problem discussed already in 1873 by Plateau In some cases their morphology and stability have 
!>■ ■ actually been elucidated experimentally in this context |3]. Other systems where minimal surfaces play a role include lipid- 
water solutions, diblock copolymers, crystallography, protein structure, or liquid crystals such as smectics |4]. They also arise as 
world-sheet instantons in string theory, for example in the semiclassical, fixed-angle high-energy limit of scattering amplitudes 

The minimal surfaces that will interest us here arise in the problem of partial wetting of a solid by a liquid 16]. In the standard 
\^ _ experimental situation, a liquid with free surface of area A (liquid-air interface) wets a flat solid plane over an area A' (liquid- 
■ solid interface). The free surface meets the solid plane along a line, called contact line, at an angle 9 which is defined locally. 
The interfacial energy is the difference E = •jA — 'j'A', where 7 is the energy per unit area (or surface tension) of the liquid-air 
-*Y> . interface, and 7' — •jsA — 7sl is the difference in surface tension between the solid-air (SA) and the solid-liquid (SL) interfaces. 
^ ' The force per unit length pushing a segment of the contact line towards the un wetted region is thus / = —7 cos 9 + j'. Requiring 
D that it vanishes gives Young's |7, 8] local equilibrium condition, 6 — aTCCos{j' / -f) . The minimal-surface problem at hand is 
i-S^ thus a problem with mixed Neumann and Dirichlet boundary conditions. In the idealized setting of an infinite liquid container 
^ and a perfectly homogeneous planar wall, there exists a simple solution to this problem: it is a planar liquid-air interface meeting 
the wall along a straight contact line. Strictly-speaking, as we will discuss in section|nj the properties of the container at infinity 
must be carefully chosen in order not to destabilize this solution. 

Two extra forces play in fact a role in the general formulation of the wetting problem. The first comes from the drop in pressure 
" ■ ' across the liquid-air interface, which adds to the Gibbs energy a volume term: £ — jA — j'A' — pV. Here p is the pressure 
difference and V the volume of the fluid. The free surfaces that minimize this energy have constant rather than vanishing mean 
curvature |8]. It is quite remarkable that the corresponding equations are (at least formally) integrable, see for instance |9|. Note 
that in the special case of an incompressible fluid, p is a Lagrange multiplier determined by the constraint that the "droplet" 
volume V be fixed. The second force that plays in general a role is gravity, which introduces an additional scale, the capillary 
length = i'^/pgY^^. Here p is the fluid mass density, and g is the gravitational acceleration. In this paper we will study 
situations where both pressure and gravity can be ignored. This is usually valid if one concentrates on length scales <C k^^, 
and considers a fluid connected to an infinite reservoir so that effectively p ~ 0. Note that the capillary length is typically of the 
order of a few millimeters, but it can be made much larger in free-fall (e.g. space-based) experiments, or if one replaces the air 
by a second non-mixing fluid of roughly equal mass density. Thus setting k ~ p 2± is a good approximation in a wide range of 
experimentally-feasible situations, and we will do so in this paper. Technically, one can further justify that gravity be ignored at 
all scales if a condition, identified below, is satisfied. 

What is in fact more questionable is the assumption of a perfectly homogeneous wall. Indeed, in most of the experimental 
setups of wetting, roughness and impurities of the solid substrate couple directly to the position of the contact line, which may as 
a result be effectively pinned. Computing the energy of a deformed contact line is thus a question of foremost importance. For 
small deformations, as Joanny and de Gennes (JdG) have shown |10], the contact line obeys non-local linear elasticity. These 
linear equations may become unstable at wavelengths comparable to some global-geometry scale, as several earlier studies have 
established Id 1.1 . The issue of non-linear elasticity, which becomes relevant for larger deformations, has been addressed only 
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recently fl?]. It could play a role f]3\ in resolving the apparent disagreement between recent experimental measurements of 
contact-line roughness 1 141. and renormalization-group calculations near the depinning transition 1 15| or numerical simulations 
flfill that were based on the JdG linear theory 1 17]. To be sure, hysteresis and other dynamical phenomena, which have at- 
tracted much of the recent attention 1 18], may also prove important in interpreting the above experimental data. Nevertheless, 
a systematic analysis should start with a thorough understanding of the non-linear and possibly non-local effects in the simpler, 
equilibrium situation. This is the problem that we will study here. 

The area of a minimal surface bounded by a given (closed) curve is simple when expressed in conformal coordinates. Non- 
linearities arise because this choice of coordinates depends non-trivially on the boundary curve, through the conformal-gauge (or 
Virasoro) conditions. In this paper we develop systematic methods for solving the ensuing non-local and non-hnear equations, 
either in perturbation theory or numerically. We focus, in particular, on the case of a planar wall, and derive simple diagrammatic 
rules that calculate the energy of a deformed contact line to any given order in the deformation amplitude. The method can be 
extended to more complicated container geometries, but the details become more involved. As a relatively simple illustration, we 
show how to extend the rules, and calculate the JdG linear theory, in the case of two contact lines lying on parallel opposite walls. 
We also describe a novel algorithm which finds the minimal surface energy with no need for surface triangulation. Finally, we 
discuss some general properties of these perturbative expansions, which bear a fascinating similarity to problems encountered 
in perturbative string theory. We hope to return to some of these questions, as well as to the implications of our results for the 
wetting problem, in a future publication. 

The paper is organized as follows: In sectionlUlwe describe our basic model, point out the need for global tadpole cancellation, 
and discuss the relation of the mixed Neumann/Dmchlet to the pure Dirichlet problem. In section|III|we give the formal solution 
of the latter problem, for an arbitrary boundary curve, in terms of conformal coordinates. This is standard material which is 
included here for completeness. In section Hvl we specialize to the case of a planar wall, derive the corresponding non-linear 
boundary equations, and express the energy in terms of their solution. We pay particular attention to the decoupling of the 
large-volume cutoff, which as we will explain is rather subtle. In section|V]we solve the boundary equations perturbatively, and 
compute the corrections to the JdG energy, up to quartic order Section fVll describes an alternative approach, using Lagrange- 
multiplier fields and leading to a simple diagrammatic representation of the perturbative expansion. The numerical algorithm is 
presented in Section lVTII In section|Vlll]we extend this to the case of two parallel walls, and calculate the quadratic interaction 
of the contact lines. Finally, in section llXl we establish the finiteness of the perturbative expansion order by order, and point out 
some intriguing directions for future work. The Weierstrass parameterization of our fluid surfaces, a calculation confirming the 
decoupling of the large-volume cutoff are described, respectively, in appendice s IaI and IbI 



11. THE MODEL 



We consider a fluid inside a tubular container x M, where M corresponds to the height coordinate z, and fl is some (a priori 
arbitrary) connected region in the {x, y) plane, with boundary dH.. Let us for now assume that the fluid surface has no overhangs 
- it can then be parameterized by the height function z{x, y). We may express the energy functional as the following sum of 
two-dimensional bulk and boundary terms: 



S <^bulk ^" *^bnry 



dxdy (^j^l + {d^z)^ + {dyzf - pz+^pgz^^ - J^^dlj'{l)z 



(1) 



where dl is the infinitesimal length along the boundary of fl. The first term in is the fluid-air interfacial energy jA, the second 
is due to the difference in pressure between air and fluid, the third to gravity, while the last comes from the fluid-solid interface. 
For convenience, we have slightly generalized the model so that the tension of this interface may vary along the container walls, 
as can be done by design. The more general case of a 7' depending on both / and z, due for instance to the presence of impurities, 
will be discussed below. For now 7' is only a function of /. 

In the absence of gravity g — 0, the minimum of the energy f is a surface of constant mean curvature, with specified contact 
angles: 
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where V = {dx, cV) and ri is a unit vector normal to the boundary dft. These non-linear equations do not always admit a global 
solution, see e.g. fly\. A necessary (but not sufficient) condition for a solution to exist is 



g = p X Area(r2) + (f dl = 
Jon 



(3) 



FIG. 1: A fluid surface bounded by a (shaded) planar wall, touching it at the position of a (pinned) contact line. At distances much larger than 
the capillary length it is flat and perpendicular to the gravitational field (left). Enlargement for distances smaller than 1/ n (right) which 
is the range of scales studied here. The unperturbed surface is a plane, making an angle 9q with the wall. When perturbed, the conformal 
parameters (ai, (T2) approach Cartesian coordinates far from the wall, as discussed in section lTVI 

This is a condition of average-force cancellation: indeed, the left-hand side of the above equation couples linearly to the zero 
mode of z{x, y), and would lead to a runaway solution if it did not vanish, the energy being unbounded in that case. By analogy 
with string theory we may refer to this as a global tadpole cancellation condition. Note, in particular, that for a homogeneous 
wall, for which 7' is constant, one must fine tune the ratio of perimeter to area so that it equals p/7'. If the average-force 
condition is satisfied, the average height of the fluid surface becomes a free dynamical parameter of the solution, analogous to 
the string-theoretic moduli. Its role must be examined with care as it threatens a priori the stability of any perturbative expansion 
at weak disorder, and may thus lead to qualitatively new behavior. 

The emergence of condition clearly originates from the neglect of gravity. If g 7^ 0, it is easy to see that the energy is 
always bounded from below and that the fluid will tend to rise such that J dxdy z ~ Q/pg, the well known capillarity effect. 
Hence if Q is non zero, one expects that the theory studied here, obtained setting g = 0, breaks down for wavevectors q < k 
(hence especially for the zero mode). However, the interesting point, discussed below, is that if one imposes Q = then one can 
safely set 5 = and obtain a theory which is well defined at all scales. This is the theory studied here. It is illustrated in Fig. [2 

Let us consider minimizing the energy in two steps: We first solve the bulk equations keeping the contact line fixed, i.e. we 
find the surface of constant mean curvature, Zfi{x, y), such that the restriction of Zh to dft is a given function h{l). We denote 
the corresponding bulk energy (or reduced energy functional in the language of lllll ) by E[h] := Ebnikizh)- The energy of the 
equilibrium configuration is then the minimum over all contact lines of 



Thus 7' plays the role of a source, and the minimum energy is just the Legendre transform of the reduced energy functional. 
If 7' were to depend also on z, the source would be field-dependent. We will comment on the subtleties of this Legendre 
transformation between the Dirichlet and Neumann problems in the concluding section. 

Let us describe the simplest configuration studied here, which consists of a semi-infinite fluid bounded by a homogeneous 
planar wall at a; = 0. We assume from now on that p = 0, and that the container at infinity has been adjusted so that the global 
tadpole condition is satisfied. The unperturbed fluid surface is then an inclined plane, making a contact angle 9q = arccos(7Q / 7) 
with the wall, as illustrated in figure 1 . We are choosing the origin of coordinates so that the unperturbed fluid surface intersects 
the wall along z = 0, while the perturbed contact line is given by z ~ h{y). It turns out to be convenient for the following to 
define: 




(4) 




V 



(5) 



If the contact line deformations are concentrated in a finite region one expects this energy difference to also be concentrated in a 
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finite region, and the outer boundaries of the container to decouple. More generally, the simple planar model of figureQlshould 
give an adequate description of the physics if all other distance scales of the system (including the capillary length, 7/p, and all 
geometric scales) are much larger than the typical deformation wavelength. We will come back to this subtle issue later on. Note 
that we have included in the energy difference the contribution, Ehmy, of the homogeneous wall. This means that E[h] should 
start out as a quadratic functional for small h{y). 

Let us briefly mention the case of impurity disorder. In this case translation symmetry of the tube is in general broken by the 
roughness of the wall. The effect of impurities can then be modeled by a variable fluid-solid tension, and the boundary term in 
([1) becomes: 



bury 



i dl [ dCi'iiX)- (6) 
Jan Jq 



The two stage minimization can then be summarized as follows. One writes: 

7'a,z) = 7^(0 + A7'(Z,z), (7) 
where 7o(0 =7 cos Oo{l) is some average or reference value, and defines the shifted functional: 



E[h] = E[h] - E[0] - 7 i dl cos9o{l)h{l) 
Jon 



(8) 

Because of disorder the impurities generate a potential for the zero mode zq of z{x, y) and the condition (|3} cannot hold in 
general. However we can still impose this condition "on average" § dl Jq{1) — and compute the corresponding E[h]. It is 
this functional which is studied here: it obeys quasi-locality and is well defined for 17 = 0. Once E[h] is known, finding the 
(equilibrium) position of the contact line amounts to solving in the second stage of the minimization: 
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E[h]- jdl J dCAj'{l,C) 



(9) 



This can be viewed as a generalized Legendre transformation, which we will not study here. The aim of this paper being simply 
to characterize E[h] in presence of an average contact angle. We will use expressions such as pinning condition, or pinned 
configuration in the following only to denote the fixed-/i conditions. 



ni. CONFORMAL COORDINATES 



Computing the area of a minimal surface bounded by a continuous closed curve v(s) is a classical problem of applied mathe- 
matics. In this section we will explain how, in conformal gauge, it reduces to a (non-linear and non-local) equation for a function 
of one variable on the boundary. Let f{ai,a2) be an arbitrary parameterization of the surface, i.e. r = {x, y, z) is the position 
of the surface S corresponding to the values of the two (a priori arbitrarily-chosen) parameters (cri, (72). We will assume that E 
has the topology of a disk, and that the parameterization is global, i.e. that there is a one-to-one correspondence between points 
of E and points in some parameter domain 2? C M^. One should of course keep in mind that, for some boundary curves, these 
assumptions may have to be relaxed. In terms of the induced metric gab = daV ■ 9f,r , the area of E reads: 

■^^ jj dCTidcra yjdetg . (10) 

This expression is invariant under any reparametrization with non-vanishing Jacobian, i.e. ai a'i(cri,o'2) and (T2 
o'2(o'i,cr2) with det [dad'h] 7^ 0. For a surface without 'overhangs' we may use this freedom to set {<Ti,a2) = {x,y), in 
which case (I10> reduces to the expression for the area used in eq. Q. This is a useful parameterization when dxZ and dyZ 
are small, but more generally the minimization of the area in this gauge leads to non-linear partial differential equations in two 
variables, which are hard to solve. 

A more convenient choice is conformal coordinates, which are defined implicitly by the two conditions : 

di'r-d2r = and dir ■ dir — ■ 82? . (11) 

Put in words, the two vector fields tangent to the surface are orthogonal everywhere and of equal, not necessarily constant, 
length. [As the reader can easily verify, the parameterization (di, (T2) = (a;, y) is conformal only in the special case of constant 
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z.] It follows from ([HJ that gab = Sab, where = dir- difis the so-called "conformal factor". Thus in this gauge the area 
can be written as 

A^^JJdaida2{dif-dir + d2r-d2r), (12) 
and the variational equations are the Laplace equations in two dimensions : 

da (^./d^g'^'dbr) = {df + di)r = 0. (13) 

The embedding coordinates {x, y, z) are therefore harmonic functions of (cti , (72), and can be written as the real parts of analytic 
functions of the complex variable w — (cti + i<72)/2 : 

x{w,w) = 2ReX{w) , y{w,w) ^ 2ReY{w) , z{w, w) ^ 2Re Z{w) . (14) 

This property of harmonic functions is very special to two dimensions. Our problem is now to determine X, Y and Z for the 
given boundary curve v(s). 

To this end, note first that if the surface is non-singular and bounded, the functions X, Y and Z must be analytic in the interior 
of the domain V. They are furthermore related by the two conformal-gauge conditions which can be combined in the 
following equivalent form : 

(9i - id2) f- {di - 182) r = (X'f + {Y' f + {Z'f = , (15) 

where the prime denotes differentiation with respect to w. This rewriting makes manifest the residual freedom of analytic 
reparametrizations of w. Such complex-analytic changes of coordinates preserve indeed the conformal condition M5\ . and 
can be used to map the parameter domain to any convenient simply-connected region in C. Let us assume, for instance, that 
I? = {w e C, \w\ < 1} is the unit disk. We write w = pe"^, and denote by r(0) :— f{<j>,p — 1) the boundary curve 
parameterized by the special conformal coordinate cj). Note that r((/)) has a unique harmonic extension to the interior of the 
disk, and thus determines unambiguously the minimal surface. This follows easily from the fact that the analytic function X{w) 
admits a Taylor expansion 

00 

XH = ^X„w;", (16) 

n=0 

SO that its restriction to the boundary has no negative-frequency Fourier modes, when identifying w" ~ e*"^". Thus, to extend 
x{(j>) to the interior of the disk, we need only split it into positive- and negative-frequency parts, x{(j>) ~ x+ {(f)) + x^ {(f). Then 
a;+ can be extended to X{w) by the replacement e**^ w, while — x^ extends to the complex-conjugate anti-analytic 
function X{w) — X{w). If x{(j)) has a zero mode, it must be split equally between the two parts. A simple calculation leads in 
fact to the following Cauchy relation between X{w) and the boundary restriction of x: 

(17) 



1 f^'^ pi-<P 
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Similar relations hold of course between Y{w) and yicj)), and also Z{w) and It is, furthermore, easy to check that since 

X(e*'^) = a; _!_((/)), the conformal-gauge condition il5i is equivalent to 

^•^ = for all e [0, 27r] . (18) 
dcp d(p 

Let us go back now to the expression il2\ for the area. If the surface is minimal, integrating by parts and using Laplace's 
equation allows to rewrite its area as a boundary integral: 



^ /•! /'27r />2'7r 



(19) 



The integrand involves the radial derivative of r, but with the help of Cauchy's equation (pdpX = —id^X, and similarly for the 
functions Y and Z) we can convert this to an angular derivative, with the result: 

a™ ^ i r . - r_ ^ ^) . 2. f: „ |F„ P ^ ,20) 
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Here fn is the Fourier transform of the function on the circle r{(f>) = ^n^™'^- For later use, we also give two alternative 
(equivalent) expressions for the minimal area: 

1 f f df df . ^ ( 4)- (t)'\ _ I f f \f{<j)) - f {(/)') f 



4^y0V d<^ d0' V 2 ; 16^7^7^, (sini(0 -(/.')) 

The first can be obtained from eq. ilOi by Fourier transform, while the second follows by a double integration by parts and the 
fact that, thanks to the ie prescription, only the cross term in the numerator contributes. Note that for suitably smooth f{(p) these 
integrals are manifestly finite in the ~> 0' region (hence the ie can be dropped in the final expression - but not if one expands 
the square). 

We have thus succeeded to express the minimal area as an explicit (non-local, but quadratic) functional of f{(f>), so one may 
think that our problem is effectively solved. This is, however, not quite the case, because the transformation from the original 
parameter of the boundary, to the special conformal coordinate 0, depends itself non-trivially on the boundary curve. To make 
this relation explicit, let us write s — /{(j)), so that f{(f>) — v(/(0)). A straightforward calculation starting from the integral 
expression ([HJ gives 

'^r. ^ v(/(^')) (22) 



d0 ^Ti J sin2(^^^) ■ 

Plugging this in the gauge condition M%\ leads to a non-linear integral equation, that can be used (in principle) to determine /(</!)) 
for any given boundary curve v(s). This is still a non-trivial task, but we have at least reduced the minimal-surface problem to 
one involving only one unknown function of a single variable. In some cases, the problem can be simplified further by using 
the residual freedom of conformal transformations to map the unit disk to a suitably-chosen domain. Such is the case when the 
contact line lies on a plane, as we will now see. 



IV. CASE OF A PLANAR WALL 
A. The boundary equations 

In the configuration of fig. 1 the contact line is restricted to a planar wall, located at a; = 0. Assuming that it has no overhangs, 
such a contact line is naturally parameterized by the height function z = h{y). We want to adapt our previous general discussion 
to this special situation. The story is somewhat simplified by using the convenient conformal coordinates (reminiscent of the 
proper- time gauge of string theory): 

ic 

X = —icw — {cr + it) , so that x — 2KeX — ct . (23) 

Here c is a positive constant, and we have traded (ai, (J2) for the lighter notation (cr, t). In imposing (I23> we have used the 
residual freedom of conformal transformations, and the fact that X is an analytic function. Note, however, that this choice of 
gauge might be obstructed globally, as we will explain in appendix A. Since the fluid surface extends out to infinity, the new 
parameter domain is the upper-half complex plane, 2? = {w G C, Imw > 0}. Later we will consider a second wall at x = L, in 
which case T> will be the infinite strip < r < L/c. The points at infinity must actually be treated with care: the right procedure 
is to first make T) finite, by bounding the fluid with outer walls, then move these outer boundaries to infinity. 

We will be here interested in surfaces that approach asymptotically the inclined plane 

fo ~ (sin 60 T, a, — cos Oq t) . (24) 
It is therefore convenient to choose c = sin Oq, and to define the difference 

Ar = f— fo, with Af=(0,y, z). (25) 

Note that the gauge condition ( I23> ensures that the first component of Af is identically zero. Since the components of both f 
and fo are harmonic, so are those of their difference Af. We can in fact write y{w, w) = 2Re Y{w) and z{w, w) = 2Re Z{w), 
where the new analytic functions are given by 

Y = Y — w ^ and Z = Z ~ i cos 6*0 w . (26) 

Following the same logic as in section|III| we also define the restrictions of y and z to the real axis, {/(cr) = y{(j, r = 0) and 
z{a) = z{a, T = 0). The extension of these functions to the upper-half plane is uniquely determined by the property that they 
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FIG. 2: The domain C and the inward pointing normal. 



should be both bounded and harmonic. Indeed, the analytic function Y must have a Fourier-Laplace expansion involving only 
positive-frequency modes: 



Y{w)=j ^ne 



2iku 



dfc 



m= I ^{Yue^^^ +C.C.) 



(27) 



since it would otherwise diverge when r ^ oo. Thus, to extend y{a) to the upper-half plane, we must first split it into its 
positive- and negative-frequency parts, y{a) = y+ic) + y-{a), then extend y+ analytically and j/_ as its complex-conjugate 
anti-analytic function. The Cauchy integral formula relating Y{w) and y{<T) reads 



Y{w) 



2^ 



dcr 



2w — a 



(28) 



The right-hand side is analytic in the upper-half complex plane provided that y{a) vanishes at infinity. Of course a similar 
formula relates also z{a) to its analytic counterpart Z{w). 

Our problem is thus reduced to that of finding the two real functions on the real axis, z{a) and y{a). These are related by the 
pinning condition of the contact line: 



z(cr) = h(^a + y{a)j 



(29) 



Furthermore, they must obey the conformal constraint (I15> . After inserting the expressions i26\ . and using the obvious identities 
y+ic) — Y{<t/2) and 5+(tT) = Z{a'/2), this constraint reads: 



d2/+ , . . dz+ 

— h I cos 6^0 — ; — 

dcr dcr 



dcr 



dz+ 
"d7 



(30) 



The pair of coupled, non-local equations i29\ and ( I30> is in principle sufficient to determine 5(cr) and y{(j), and hence also the 
complete shape of the fluid surface. In the following sections we will discuss how to solve these equations numerically, or by a 
series expansion in powers of h{y). First, however, we must express the energy in terms of the two boundary functions z{a) and 



B. Expression for the energy 



The area of an infinite fluid surface is, clearly, infinite. However, for a localized deformation of the contact line, i.e. for 
h{y) — > when y ±oo, we expect the difference in area, Amm = Amin[h] — -4,nin[0], to be finite. To calculate this 
difference, we will introduce as a physical cutoff a tubular container C = O x M, with fl a. rectangle of size x Ly in the {x, y) 
plane. We define the associated characteristic function 



Qc{r) :=e(x)eU + 



1 if r G C 
otherwise , 



with Q{a) the usual Heaviside step function. The difference of the areas then reads: 



gab 



(31) 



(32) 
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where, after evaluating the right-hand side, we should take the limit L^, Ly oo. Note that cutting off directly the parameter 
range could give a wrong answer, because the same value of (ct, t) need not correspond to the same value of (x, y) on the planar 
and on the deformed surface. 



Expanding the integrand of eq. (I32> in powers of Ar, and using the fact that 

ec(ro + Af) = ec(n)) + 5c{rQ) n ■ Ar + ■ 



(33) 



where 6c is the delta function localized on the boundary of C and h is the inward-pointing normal unit vector, leads to the 
following expression for the area difference: 



{daAf-dbAr + 2daAr- dbfo) S"" + \h ■ d^ro\-' n ■ A 



V 



(34) 



dV 



Here V = [0,Lx/ sin^o] x [~^v/'^t Ly/2] is the parameter domain defined by the condition Scl'ti) = 1> and dV is its 
boundary. The last term in the above equation accounts for the fact that the cutoff corresponds to a container in physical space, 
rather than in the space of parameters (cr, r). The factor \n-d±fo\~^, with d± :— h-d a derivative in the direction normal to dD, 
is the Jacobian that arises upon converting Sc{ro) to a (5-function in parameter space. The neglected terms involve higher powers 
of fi ■ Ar, and one or more partial derivatives. They vanish on the outer boundary, provided A?- ^ at infinity, and on the x — 
wall where n ■ Ar = for our choice of gauge. Note that in deriving expression i34l we used the equality daro • dbro (5°'' = 2, 
which follows easily from ( I24t . 

Using Stoke's theorem and Laplace's equation we can express all the terms in i34i as boundary integrals. 



-An 



dD 



-Ar • d±Af + Ar - d±rQ — \n ■ dj_rQ\ ^ n ■ Ar 



(35) 



Let us consider first the a = ±Ly/2 boundaries: since h = =p(0, 1, 0) there, the last two terms cancel exactly one another, while 
the term quadratic in Ar does not contribute as long as Ar ^ at infinity. This term does not contribute, for the same reason, 
at the T = Lx/ s'mOo boundary. Finally, at both r = and t — L^/ s'mOo we have n ■ Ar = 0, since n = ±(1, 0, 0) on these 
boundaries and, with our choice of gauge, At^^ (0, y, z). Putting all these facts together we obtain: 



Amin = J da ^--{ydr 



y + zdrz) 



+ cos 9o z 



cos Uq z 



9o i" 



(36) 



where the limit Ly ^ oo has already been taken on the right-hand-side. The r-derivatives in the first term can be converted to 
cr-derivatives with the help of the Cauchy equation. As for the last two linear terms, they cancel because z is harmonic (both are 
proportional to the same fc = Fourier mode). Thus the difference of the areas reads: 



da dcr 



+ c.c. 



(37) 



Although this calculation is correct, the cancellation of the linear terms is, from the physical point of view, rather misleading. 
It involves two opposite walls which are infinitely far apart in the ^ 00 limit, and looks therefore highly non-local. A 
physically more significant cancellation occurs in the energy functional E[h], which (as explained in section |ll} receives a 
contribution from the fluid-solid interface: 



-'bnry 



dV 



d/7'z = —7 COS 6*0 / da ^z{l + daij) 



=L^/ sin 00 



(38) 



The second equality can be understood as follows: the unperturbed planar surface meets the x = 0, x = Lx and y = ±Ly/2 
walls at angles equal to , vr — 6*0 and tt/2, respectively. Young's equilibrium condition thus requires that, in the absence of 
impurities: 



7- 




7 cos 9q 



for X — , 
for y = ±Ly/2 , 
for X — Lx ■ 



(39) 



Furthermore, along the first and the last wall the invariant length is dl — dy — [1 + day) da. Dropping the quadratic term at 
X = Lx, since both y and z must tend there to zero, gives the advertised equation ( I38t . Adding this to 7-4inin leads to our final 
expression for the energy: 
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Note that the linear terms cancel here separately on each wall, and that all the contributions to the energy are "quasi-local". Thus 
the large-volume cutoff decouples, as expected, in the calculation of the energy (but not of the separate contributions jAmin 
and f bnry)- The only restriction on the cutoff is that it should not destabilize the unperturbed planar surface. We confirm these 
claims by a calculation in appendix B, which includes as an extra control parameter the inclination angle of the outer wall. 

For later use, we will also need the expression of the energy in terms of the Fourier components of y{a) and z{(t). Using 
eq. ( I27t . and doing some straightforward algebra leads to 

f°° dk / - ~ ~ \ 

E[h]^jJ — ki^\Yk+ icos9oZk\^ + sin" eo\Zk\^) . (41) 

Note that the energy is quadratic in y and z, where the function y{a) = a + y{cr) relates the natural parameterization of the 
contact line, to the conformal parameterization in terms of a. As was explained in the previous section, the problem is non-linear 
because this change of coordinate depends explicitly on the pinning profile. 



V. PERTURBATIVE EXPANSION 



The pair of equations ( I29l l and OOt cannot be solved, in general, in closed form. However, if the contact line is deformed only 
"slightly" (this will be made more precise later), then y and z should both be small. We may therefore expand the right-hand 
side of eq. J29t in a Taylor series. 



n=0 



d"fe(a) yjaf 
dcr" n\ 



(42) 



where both y and the derivatives of h are now evaluated at the argument a. Furthermore, solving the quadratic equation ( I30> for 
dy+ /da and integrating gives: 



da' 



4 VdCT 



i cos 9o 



da' 



(43) 



Note that we have picked the solution of the quadratic equation that vanishes for z+ 0, and we have also fixed arbitrarily the 
irrelevant (complex) integration constant. Since z is small, we may expand the integrand on the right-hand side to find 



y+{<^) 



Tn — 1 



da' 



[(-1) • 1 •3-5---(2n-3)] X 

dz+ . r, „ /dz- 



da' 



— oo 
2 



/dii \2 . dz+ 



(44) 



Equations ( I42t and ( I44t can now be solved iteratively as follows: one starts with the lowest-order solution of the first equation, 
z{a) = h{a), and inserts it into the second one to find j/+ — —i cos 6*0 Inserting the result in eq. Hli gives z at quadratic 
order in h, and from eq. ( I44t we can obtain y to the same order Iterating the procedure gives, in principle, the solution to any 
desired order in the pinning profile h. 

In order to write the answer in a compact form, we introduce the following notation. If f± {a) are the positive- and negative- 
frequency parts of any real function f{a), then 



f = f+ + ,U and i f ■.= !+- f- 



(45) 



where the second equality defines the dual function f{a). Note that / and / are both real - this follows from the fact that 
/_ = (/+)*. Now the first few orders in the expansion of the solution read: 



z — 11 + COS Oo — n 
da 



cos^ 00 d'^h 
~ d^ 



dh 
d^ 



sin^ 9o 



dh^ 



= —i cosOq h-^- — sin 



dh+ 



da' — i cos^ 00 { 



oo 

dh 
d^ 



da' 



da' + i cos^ 6*0 

,3\ 



dh. 



\da 



+ c.c. 



0{h^), (46) 



(47) 



where we have stopped at one order lower in the expansion of y for a reason that will become apparent in a minute. It will be 
useful to have also at hand the Fourier transforms of these expressions. Noting that 



1 — ^ k 

f+g+ + f-9- = 2^f9 -fa) and i = fk r^^ 



(48) 
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we find after some straightforward manipulations: 

iYk = hk cos 00 - 
Zk = hk+ / hk^hk^ kik2 



hkihk2 kik2 
cos 6*0 



sin^ ^0 



e(fcifc2) 



cos^ 6o 



\k2 



hk^hk^hk^ kik2k3 



siu^do , , cos^ 6*0^1 cos^ 00 (^2 + fcs) 

■fc)(fc2K3) 



k2 + ks 



2|fc2fc3| 



1^3 life 



(49) 



(50) 



Here the integrals run over all kj, with normalization / dfcj/(27r) and the condition that ^ kj — k. The step functions Q(kikj) 
force the two momenta to have the same sign, and we have assumed that k is positive. Recall that Yk enters into the expression 
( 14 H for the energy through the combination 



i{Yk +« cos 6*0 Zk) 



siv? 00 



hk^hk2 fcifc2e(fcifc2) + 0{h^) 



(51) 



Since this starts out quadratically in h, the cubic corrections contribute to the energy at 0{h^). This explains why we have 
truncated the expansion of y at one order lower than the expansion of z. 

Inserting i5Q\ and (15 1> in (I41> . and doing some straightforward manipulations, leads to the following expression for the energy 
of the deformed contact line at quartic order: 



E[h] ^E2+E3 + E4 + 0{h^ 



where 



£'2=7 sin Oq 



dk 
2tt 



(52) 



£^3 = 7 cos 6*0 sin 9o / hk^hk^hk 



\k1\k2kz 



-7 cos 00 sm 



% J hk. 



(53) 



hk^hk^hk^hki kik2kzk4. 



. 4. re(fcifc2)e(A:3fc4) 

sm 00 \ - 



\ki 



2ki eikski) 
N {ks + ki) 



■4)/ 



sm 00 cos 00 



kiki 



k2 



k^ 



Ifefcsfel |^i||fe||fci + fel 



(54) 



The integrals in ( I53> and ( I54t run over all kj with the condition that ^ kj = 0. As a check, note that for Oq = tt/2 the energy is 
invariant under reflection, h ^ —h, of the contact line. Note also that the expressions multiplying J| hkj inside the integrals are 
invariant under the combination of complex conjugation and change of si gn o f all the momenta, consistently with the fact that 
E[h] should be real. The expression for E3 agrees with the one derived in Il3l by a different method. 

The Joanny-de Gennes linear theory 1 10] corresponds to the leading term of the above expansion. Comparing E2 with the 
energy of an elastic rod, _E ~ J k'^\hk\^, one notes a softening of short-distance modes, and corresponding hardening of long- 
distance modes, due to the interactions mediated by the surface. In real space, the JdG energy can be written as (see the discussion 
in sectioninU: 



Eo 



47r 



dada 



{a- a' + leY 



(55) 



This quadratic, non-local functional has appeared in a variety of other contexts, e.g. in simple models of quantum-mechanical 
dissipation |20, 21]. Note that E2 is invariant under SL(2,M) transformations, i.e. under conformal transformations that preserve 



the upper-half complex plane, a 



j,h h, with a, b, c, d real and ad — bc^ 1. The full energy is not only translationaUy- 



invariant, but it also transforms covariantly under rescalings of the physical space: 



E[h'^^^] = E[h] if h'^^\y) = Xh{X-^y) 



(56) 



This implies that the perturbative expansion is really an expansion in derivatives, as should be expected from the fact that the 
classical problem has no intrinsic length scale. We will return to this point later on. 
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It will be useful, for comparison with the following section, to rewrite the quartic contributions to the energy differently. 
First, we note that the two terms multiplying sin"* 9q are equal up to a factor of —2. This follows from the following chain of 
replacements, which are allowed upon symmetrization of the integrand: 



2ki Qihki) 



<Si + S2) 



(fci + k 



\ki + k2\ 



Here Sj = kj/\kj \ is the sign of the momentum kj, and in the second step we have used the fact that the sign of (fci + fc2) is the 
same as the sign of either fci or k2, since the expression is multiplied by (1 + S1S2) = 28(fcifc2)- Likewise, one can justify the 
following replacement: 



|fci| (fca + fc4) 



kik^ 



|fcl||fc4||fcl +fc2 



and 



^2 



fcl^ 



|fcl||fc4||fcl +fc2| 



(k^ + k^) — {ki + k^) 

2|fci||fc4||fcl +fc2| 



flitting all these facts together, using that kj — and doing some straightforward rearrangements leads to the following 
alternative expression for the quartic energy: 



e(fcifc2)e(fc3fc4) 

|fcl+fc2| 



cos^ 6, 



kiki 



fc2fc3 



fc2fc3fc4| 



|fcl||fc4||fcl 



(57) 



This somewhat more economical expression will be easier to compare with the diagrammatic expansion, to which we will now 
turn our attention. Note that the expression for E4 in the particular case Oq — 7r/2 was also found in L13J using the perturbative 
solution of the non-linear equation (not using conformal coordinates). It is possible, though cumbersome, to extend the method 
to ai-biti-ary 6*0 10 . 



VI. DIAGRAMMATIC METHOD 



The perturbative expansion of the energy can be organized efficiently by using a Lagrange-multiplier field to impose the 
pinning constraint of the contact line. One starts with the following variational principle for the area: 



Anin = cxtr^(a,r) , with A{a,r) = d^a^/detg — / dsQ;(s) z{s) — /i(?/(s)) 

JJv ' Jdv L 



(58) 



Here s parameterizes the boundary of the domain V, and a is a Lagrange-multiplier field that transforms under reparametriza- 
tions such that a{s) ds remains unchanged. Since A{a, f) is reparametrization-invariant, we are free to choose the conformal 
gauge and to set x = sin 6*0 t as before. Thus V is the upper-half plane r > 0, and we may choose s ^ a for the boundary 
parameter We also define y — a + y and z = — cos 9qt + z , and we subtract from A the area of the flat fluid surface. This 
gives ^inin — extr A, where 



A{a,y,z) = ^ II (dayd'^y + dazd^z) - j a {z - h{a + y)) - coaOo z 



(59) 



The last term in the above expression comes from the cross term daZo = — cos 6*0 drZ in the area difference. This is a 
total derivative, which is why it only contributes a boundary term. Note that, in the light of our discussion in section IIVI all 
contributions from the boundaries at infinity have been dropped. This is legitimate since we are ultimately interested in the 
energy rather than in the area of the fluid surface. Alternatively, one can view A{a, f) as an action and consider the path 
integral over the fields r and a f2§\. Since we are doing only a tree level calculation there is no need to worry about Fadeev- 
Popov ghosts, which would be important for the study of thermal or quantum fluctuations. Fluctuating surfaces |23] are beyond 
the scope of the present study. 

It looks, at first sight, rather odd that in the above formulation the conformal-gauge conditions are not explicitly imposed. The 
extrema of A{a, y, z) should therefore obey these conditions automatically. To see why, note that the variation of ( I59I I leads to 
the boundary equations: 

drjj = a{a) h' {(J + y) and drZ = —a{(T) + cos 60 at t = 0. (60) 
From the above boundary equations, and from the pinning constraint z = h{y), we deduce: 

d^-f = (sin 6*0 , a{a) h'{y) , -a{a)) and d^r = (O , d^y , h'{y) d^y) . (61) 
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Thus, on the boundary, the condition dar- drr = holds. This implies that the function d^r- d^r, which is analytic in the upper 
half plane and vanishes at infinity, has zero imaginary part on the real axis. From the Cauchy-Poisson integral formula |24] we 
conclude that it vanishes everywhere, so that the conformal gauge conditions M5\ are indeed satisfied. 

In order to develop simple diagrammatic rules, we first solve the harmonic equation for the "bulk" fields keeping their restric- 
tions to the boundary, y{a) :— y{a, 0) and z{a) :— z{a, 0), fixed. As has been akeady discussed, this leads to the replacement 



,~ , Ay+ 1 



do- 



da 



\k\yk y~ 



(62) 



and likewise for the field z. Next, we solve the linear equations for z{a), thus eliminating it entirely from the expression (15 9> . 
The new variational functional, expressed in terms of Fourier components, reads 



A[a,y) 



\k\yk y- 



Ifc 



a-hHk 



(63) 



where ak = cos 2Trd{k) + 5ak , and Hk is the Fourier transform of H{a) — h{a + y{(j))- This result also follows if one uses 
the path integral formulation and integrates over the fields y and z in the bulk. More explicitly 



Hk = hk + j iki hki yk2 + ^ Jiiki fhk^ yk^ Vki + 



(64) 



where the integrals run over ^ kj = k. The extremum of the functional ( I63t can be computed by summing tree-level diagrams 
of a 1 -dimensional field theory. The 1 -point function and propagators read: 



C := {ak) = cos6'o27r(5(fc) 

:= {6ak5a-k) = \k\ 

^\f\J\j {yky-k) = -TTT 



(65) 



while the first few vertices, deriving from the last term of ( I63t . are as follows: 

• := hka-k , S ■= /ifci yfe2 "-fei-fc2 I 




hki yk2 Vks OL-ki-k2-k3 



(66) 



Note that all of these vertices are proportional to the amplitude of the pinning profile. Furthermore wiggly lines, corresponding 
to the field y, can only terminate on another vertex in a vacuum tree diagram. Thus only a finite number of vertices contributes 
to a given order in the expansion in h. Solid lines corresponding to the Lagrange-multiplier field a may end at the tadpole 
{ak) — cos ^0 27r5(fc), which carries no extra power of h. Note also that at the vertices momentum is injected by hk, which has 
to be taken into account for momentum conservation. 

Using the above diagrammatic rules, one can compute any desired order in the expansion of E[h]. This is obtained by 
multiplying the extremum of ( I63> with 7, and then subtracting the hnear contribution of the wall, f bmy — l' J h — j cos Oq Hq 
(see section HVt . This contribution cancels precisely the tadpole diagram 

• C = cos 9o ho 



(67) 

in agreement with the fact that the unperturbed, planar fluid surface should be stable. Denoting by En the nth-order term in the 
perturbative expansion of E[h] one finds: 



Order 2: 



so that 



\k\ hk h^k , 



1 f P 

-COS^^o / jT-rhkh^k 

2 Jk \k\ 



(68) 



E2 = 1 



= ^ sin^ 6*0 / |fc| hkh^k 



(69) 



which is precisely the Joanny-de Gennes quadratic energy. 
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Order 3: 



cosfo / — 71—, — hk.hk,hk3 



(70) 



1 



(71) 



These two contributions are of the same form. To see why, one must symmetrize the integrands over all permutations of [123], 
and then use the identities that follow from momentum conservation: 

— S1S3 fc2^ + pcrnis = Si S3 (fc]^^ + fcifcs) + perms = — S1S3 /ci/c2 + perms , 

where Sj — kj/\kj \ is the sign of kj. Thus, the sum of the two diagrams gives 

S3 = 7 



which agrees with the calculation i53\ of the previous section. 
Order 4: 



7cosfc'o sm 6^0 / — rr~\ — 'ifei'ifc2'ifc3 



(72) 



1 f \k1\k2k3\k . 

2 J Ifci + k2 



(73) 



= — cos Wn 
2 ' 



fci fc| fcs I A:4 1 
|fci||A:3| 



(74) 



2 1 fcifc|fc3|A:4| , , , , 

cos 6^0 / -n— rn — —r^ hk^hu^^hu^hu^ 



|fcl||fcl+fc2l 



(75) 



f\J'\J\3( 



= — cos 
2 



NN 



(76) 




= -cos t^o 
6 



kik2k-ikl , , , , 
fci fc2 K3 



(77) 



k^kokrik^. Till 



(78) 



Note that the power of cos Oq corresponds to the number of "hooks" of a given diagram. For cos6'o — only the first of these 
diagrams contributes. Using the replacements 



S1S4 



(Sl +S2)(S3 +S4) _ e(fcifc2)e(fc3fc4) 



Ifci + fco 



4|fci + fc2 



Ifci +fc2| 



one can check that this diagram agrees with our previous result ( I57l i. To show that the expression also agrees with E4 = 
7 [(EH + GH) + (ESj + + (|77j + ijTll)] for arbitary contact angle 6*0 we proceed as follows: first diagrams (|74|i and (|77|i 
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can be combined to reproduce the second term in (15 7> . To see why, one must replace ^k^ = — ■|fc4(fci +^2 + ^3) 
—kik4 in the integrand of the diagram ( I77> . Secondly, one can show that for cos^o — 1 the sum of the remaining diagrams, 
IISJ + IIB + m + EH), is exactly zero. Indeed, writing the integrands of these diagrams in the order of their appearance we 
find: 



kik2k^ki 



\ki \ 1^4] \ki + k2 

kik2ky,ki 



[kiki + 2^:2^:4 + {ki + k2f + /S2fc3] = 



kik2k3k4 



|fci| Ik^l \ki+k2 



\ki\ \ki\ \ki + k2 



{k2ki - kiks) 



kik2k3k4 



\kl\ |fc4| |fc3 + fc4 



(fc2fc4 - kika) , 



[kiki + 2k2ki — (fci + fc2)(fc3 + ki) + A;2fc3] 

(79) 



where in the first and the third equality we have used the conservation of momentum. Since the last two expressions are equal 
they can be replaced by their average. The result is antisymmetric under the exchange of 1 with 4 and 2 with 3, so after 
multiplication with hk^ hk^ hk^ hk^ it gives zero as claimed. We are thus free to subtract this vanishing expression times ^ cos^ Oq 
from the sum of all diagrams. This removes the contributions il5\ and j76> . and changes the coefficients of j73> and ( I78> to 
those of the corresponding terms in i5H . This completes the proof that equation (|^} agrees with the diagrammatic calculation 
of the energy. 

The diagrammatic expansion can be extended to higher orders. As an illustration, let us consider the case of a perpendicular 
contact angle, in which case the tadpole vanishes. The symmetry under h —h guarantees that only even powers appear in the 
expansion of E[h]. The sixth and eighth order terms are given by: 



= 1 



\ki + k2\\k^ + ke 



(80) 




where 



Ifcl|fc2fcg|fc4|fc||fc6|fc7|fc8| 
l^'l + fePl + + ^3 + ^411^7 + fcs 



(81) 



and 




|fcl| fc2 |fc3| kj [fcsl fcg \kj\ kl 
+A:2||fc3 + fc4||fc5 + fc6| 



(82) 



We will comment further on these results in the final section. 



VII. NUMERICAL ALGORITHM 



As was shown in section Hvl the problem of determining the deformed fluid surface with a pinned contact line on a planar 
wall reduces to that of solving the pair of equations for the real functions y{cr) and z{a): 

~zia)^h(a + yia)) and ^ + -os^o ^ = - (^J - (^J • (83) 

We recall that f± are the projections onto positive (negative) momentum Fourier components of the function /. Continuing /+ 
as an analytic function, F, in the upper-half w — [a + it)/ 2 plane, determines the unique harmonic extension of the function, 
f — 2 KeF. The conformally-parameterized minimal surface is 



{x, y, z) = (sin 6*0 t, u + — cos Oqt + z) , 



(84) 
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FIG. 3: Convergence of the energy En for a Gaussian with almost maximal amplitude, as function of iteration n. Also plotted are the 
perturbative results E2 = 0.00951444, and E2 + E4 = 0.00837429. The second plot shows convergence on a logj^Q scale. Convergence 
improves considerably for smaller amplitude of h. 
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FIG. 4: A periodically repeated Gaussian for h{y), with no 0-mode: f^^ dyh{y) — 0. The corresponding functions y{a) and z{a) on the 
boundary are also given. One remarks that z almost has a cusp-like singularity at cr — 1/2. Further increasing the amplitude of h, z will 
develop this cusp, which signals the breakdown of our parameterization. 



and it has a total energy given by eqs. ( I40l l or J41t . 

The equations ( I83> can be solved by iteration, starting with the initial configuration 

y(o)('^)=0 , 7S''\(j) ^ h{a) . (85) 

Let y^"^ and z^"^ be the solution of the equations after n steps. We extract z^J^'' by doing a double Fourier transform. Plugging 



the result in eqs. ( I83> then gives the improved values of the unknown functions: 



~(n+l) 



da' 



dz 



(n) 



da' 



— I cos Wo - 



da' 



1/2 



(86) 



Using J41> yields an approximation iJ„ to the true energy E^. We have used this iterative algorithm for h,(y) = e f{y) with 
f{y) various trial pinning profiles. We found that it converges rapidly to the perturbative result for small e, and that it breaks 
down at some critical e where the function y{a) stops being monotonic. We believe this signals a coordinate, rather than a real 
geometric singularity, as is observed in section|A| If so, it would be very interesting to develop alternative algorithms that could 
circumvent this problem. 

On figure|3]we show the convergence of the algorithm at 6* = n/2 for a profile h{y) given on figure|3 together with the 
corresponding functions y(r) and z{t). One sees on figure |3 already the emergence of a linear cusp at the tip of z{t = 1/2), 
which signals for larger e the break-down of the algorithm. 



VIII. INTERACTION BETWEEN CONTACT LINES 



As another application of the general approach, we will now calculate the interaction between the two contact lines of a liquid 
surface bounded by parallel walls. For an analogous calculation in open string theory see reference 125.1 . Suppose that wall 1 is 
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located at x = 0, wall 2 at x = L, and let = —72 = 7 cos . In the absence of impurities the equilibrium configuration is 
thus an inclined planar surface making a contact angle 6*0 (respectively tt — 6*0) with the first (second) wall. We use conformal 
coordinates and set x = sin Oq t, so that the parameter domain is the infinite strip < r < L/ sin 6*0 = tq. Repeating the same 
steps as in the previous section leads to the following variational functional for the minimal area: 



0<T<To 



{davd^-y + dazd°-z) - / ai {z - hi{a + y)) - cosOaz 



T=0 



a2 {z ~ h2{(T + y)) + cos^o z 



(87) 



Here hj{y) (for J — 1, 2) are the deformations of the two contact lines away from their equilibrium configuration, and aj are 
the corresponding Lagrange-multiplier fields. The minimal area difference is ^min = cxtr^ , where one must extremize A 
over the bulk fields y{a, r) and z(cr, r) and the boundary fields aj{a). 

First we solve the harmonic equations for y and z, keeping their values on the boundary fixed. Let, for example, y{a, 0) = 
yi{<j) and y{a, tq) — §2(0') ■ Eliminating the field in the interior gives 



day d^y 



0<r<ro 



An ^ 

J, J' 



(88) 



where 



-\ogsmh'[^{a ~ a')] if J = J' , 



-logcosh2[^(a-a')] if J /J' 



(89) 



One way of establishing this formula, is to start from the analogous expression for the unit disk, eq. (I2H . and then apply the 
conformal transformation that maps the unit disk onto the infinite strip: 



V = (T + IT 



TO, /.I 
— log I — 
TT V 1 



w = pe 



tanh 



2^ 



(90) 



Notice that the two unit-radius semi-circles, p = 1 and e [0, tt] or </) G [tt, 2tt\, are indeed mapped onto the two boundaries of 
the strip, Imw = or Imu = tq. On these boundaries 



log sin 



= logsinh 



7;—iv-v) 

ZTq 



log cosh ( — ) — log cosh ( 



To / 



To 



(91) 



up to an irrelevant constant. The terms depending only on v, or only on v', will drop out when inserted in the double integral 
nil . Setting finally v — v' — a — a' (or v — v' = a — a' — ztq) for points on the same (or opposite) boundaries of the infinite 
strip, leads to the expressions ( I88> and (I89> . as claimed. An alternative derivation of this result using the massless propagator on 
the strip is 
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<T|/ro^ (93) 



These formulae agree with ( I89l l up to an irrelevant constant. 

It will be useful to write these expressions in Fourier space. This can be done by using the identities 

^ (-1)" TT 



1 TT 

= — coth(7r6) , and > , „ t — ,.,/, \ 

n— — 00 



(94) 



To lighten the notation, we will suppress the label of the boundaries, and use boldface letters for the corresponding vectors and 
matrices. Thus y will stand for the two-component vector (yi, ^2); and G for the 2x2 matrix-valued kernel Gjjt. With the 
help of the above formulae one finds: 
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dy 
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where t indicates the transpose of a vector, and 



G(fc) 



k coth(rofc) 
-fc/sinh(Tofc) 



— fc/sinh(Tofc) 
fccoth(Tofc) 



(96) 



Since det G(fc) = k^, the inverse matrix takes also a simple form: 



coth(ToA;) 
l/sinh(Tofc) 



As a check note that in the limit of an infinitely-wide strip (L tq 
kernel for two separate, half-infinite planes. 

Returning to the variational functional ( I87> . it can be replaced by 



l/sinh(Tofc) \ 

coth(rofc) / 
► oo) one finds G{k) ~ 



(97) 



This is indeed the 



y^G{k) y_fc + i^'Gik) i-u] + (a,. • H u - ■ i-k) 



(98) 



where a = (ai, 0:2) is the vector of Lagrange-multiplier fields, (5a = a — (a) ~ (ai —cos 0^2 +cos ^o)^ and Hfc is the Fourier 
transform of the (vector of) composite fields hj{a + y,j{a)). Solving the linear equations for z, and inserting the solution in the 
above functional gives 



A{a,y) = ]^ jjkG{k)y^k-\ jjakG{k)-Hs.^k + ^ a,. • H^^ . 



(99) 



We can now read off the Feynman rules that generalize the ones of the previous section. The propagators and 1-point functions 
for the vector fields are: 



C := (afc) = (cos6'o,-cos6'o)^(fc) 

:= (5afe(5al,,) = G(fc) 



(100) 



The vertices do not mix fields on opposite boundaries, and are thus two copies of the vertices in ( I66> . 

Using these rules we may calculate the energy to any desired order in the /i-expansion. The leading, quadratic energy that 
generalizes the JdG result reads: 
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(101) 



Since both terms inside the integral are positive-definite, it is energetically favorable for hi k and h2.k to have the same phase. 
Thus the interaction between the two contact lines is attractive. Note that if we fix hi and allow /i2 to freely adjust, we find that 
the minimum of the energy is obtained for 



h2(k) 



hi(k) 



cosh(fc_L/ sin ^o) 



The energy for given hi and free /i2 thus reads 
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(102) 



(103) 



In the limit of L ^ cxi, we recover our previous expression i52\ as expected. 

Taking the same limit in dlOU shows that the interaction decays exponentially, as ~ exp(— 2fcL/ sin6'o). This exponential 
decay also applies for fixed L and very small contact angle, since the actual separation of the (unperturbed) contact lines is 
L/ sin 00- In the opposite limit of a thin strip, or equivalently of very long-wavelength deformations, we find: 
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The leading term has a simple geometric interpretation: It is proportional to the increase in area of a planar strip, whose bound- 
aries undergo a relative displacement hi — h2 along the walls, with which it made initially an angle 6q. For hi — h2, the next 
term in the above quadratic energy corresponds to an elastic rod with effective tension ^^s = iL sin ^o- This has also a simple 
geometric interpretation: The rod is in fact a thin strip, of width L / sin ^o, which is deformed by an amount hi{y) sin 6*0 in the 
transverse direction. 



IX. DISCUSSION 



In the previous sections we have shown how to calculate the energy of a deformed, almost rectilinear, contact line to any 
desired order in perturbation theory. We would now like to discuss some general properties of this expansion. One important 
point is that perturbation theory is quasi-local, i.e. the total energy is concentrated in a region of size equal to the typical 
wavelength of the deformation. We indeed saw that, as long as the large-volume cutoff has been fine-tuned so as to cancel the 
global tadpole, it decouples from any localized perturbation. One would expect the same to be true for all other geometric length 
scales of the system, such as the wall's inverse curvature. If this is true, at sufficiently short distances perturbation theory should 
be scale-covariant, as was pointed out in section|V] In momentum space, the scaling symmetry ( I56> reads: 

£;[/i(^)] = E[h] for h^^^ = hxk ■ (105) 

Inspection of eqs. (I68ll-(I82|| shows that this indeed holds at each order of the expansion, and even for each individual diagram. 
Note, in passing, that the scaling symmetry does not imply conformal invariance, as would have been the case if the one- 
dimensional theory were truly local. 

Finiteness of the JdG quadratic energy requires that 

khk , for both k ^ oo and fc ^ . (106) 

In other words, h{y) must be continuous everywhere and finite, and it must vanish as y — > ±00. A more stringent condition is, 
in fact, required to prove ultraviolet finiteness at all higher orders. It reads 

k^hk -^0 for fc ^ 00 h[^'^ ^0 for A ^ 00 . (107) 

Stated differently, the profile function h{y) must also have a continuous first derivative. That this is indeed necessary follows 
by considering for instance the "comb" diagrams, the first few being ( I78t . ( I80> and ( I81> . As the reader can check, a power 
fall-off slower than would make the comb diagrams with a large enough number of vertices diverge. To show that this 
condition is also sufficient, it is convenient to assign the scaling dimensions [k] ~ 1 and [hk] = — 2 to the factors entering 
in a diagram. Because k^hk — > at high momentum, the degree of divergence of any partial integration is bounded from 
above by the corresponding scaling dimension, in which one only counts elements that depend on the integrated momenta. The 
scaling symmetry ( I105> implies that the overall scaling dimension of any tree diagram is —2, so there is no divergence from the 
integration region where all the momenta go to infinity. Keeping one (or more) of the momenta fixed amounts to removing from 
the counting a factor dfc /c™ hk, and at most one solid and m curly propagators that emanate from the corresponding vertex. This 
can only lower the scaling dimension, so all the partial integrations are also ultraviolet finite, q.e.d. 

Infrared finiteness is trickier to establish diagrammatically. Condition ( IIO6I 1 suffices to ensure that there is no divergence when 
the momenta flowing into individual vertices go to zero. The dangerous diagrams are, however, those for which such momenta 
add up to zero along some curly line. Inspection of the expression i54\ shows, nevertheless, that the result is finite up to quartic 
order, thanks to the Heaviside functions that multiply such dangerous terms. To prove finiteness at all higher orders, it is more 
convenient to go back to the pair of classical equations ( I29> and ( I30l l. Let y'") (a) and z*^"' (a) be the solutions of these equations 
at nth order. It is then straightforward to check that, if these functions vanish at y ^ ±00 for all n < N, they will continue to 
do so for n = N + 1. This is in turn sufficient to guarantee the infrared finiteness of the energy at all orders. 

What about non-perturbative effects? To fix ideas, let h{y) — e f{y) with f{y) a given profile function, and e the parameter 
controlling the perturbative expansion. One expects that the radius of convergence of this expansion is finite, since at large 
enough e the solution to eqs. M9\ and ( I30> should stop being analytic. This could signal either one of the following two things: 
(i) that our parameterization is singular, or (ii) that the surface develops real geometric singularities or that there is a change 
in topology. It would be very interesting to find some general criteria which could distinguish between these two possibilities. 
Note that a topological transition may occur if it is energetically favorable to drill two holes in the fluid surface, and to replace 
the corresponding disks by a cylinder. In any case, the following simple (though rather crude) linear bound 



E[h] <j \h\+y h 



y 



(108) 



19 



shows that the energy of a pinned contact line stays finite. Furthermore, localized microscopic perturbations always have a 
vanishingly-small energy, and should decouple from the physics at longer scales. 

This brings us to our final remark ll26ll : as was explained in section |lll the purely-Dirichlet minimal surface problem is 
related to the mixed Dirichlet-Neumann problem, relevant for capillary phenomena, by a Legendre transformation. A Legendre 
transformation looks at first sight rather benign, but it drastically modifies the nature of perturbation theory. This is best illustrated 
by the following spectacular phenomenon |27|: A wedge in the tubular contour dil of section|lll with opening angle less than 
\tt — 29q\, is a local geometric obstruction which forces the capillary surface to develop a second sheet. This has been observed 
in micro-gravity experiments. Notice that the wedge can be of microscopic transverse size, but it should extend to all values of 
the height coordinate z. It is the latter assumption which is responsible for the apparent non-decoupling of short-distance scales. 
The story is reminiscent of the role played by wormholes in theories of quantum gravity. This analogy, as well as the possible 
impact of wedge defects on the problem of wetting, deserve further investigation. 
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APPENDIX A: WEIERSTRASS PARAMETERIZATION 



In section llVl we have parameterized the minimal surfaces in terms of two functions Y{w) and Z{w), which are related by the 
conformal-gauge condition ( I30> . The parameterization is global provided the two functions are analytic everywhere in the upper- 
half complex plane. This is indeed the case in perturbation theory, but more generally, for a given analytic Z, the solution of 
i30\ will not give an analytic function Y. A constructive solution of the conformal-gauge condition, that guarantees analyticity, 
is given by the Weierstrass representation 



(w) = dvf X \ - 




(Al) 



where f{v) and g{v) are holomorphic functions in the upper-half plane. To go to the special gauge eq. i23\ of section Hvl one 
sets 2/ = ic/g. The surfaces are then parameterized by a single function: 




(A2) 



Clearly, this special parameterization is non-singular if and only if g has no zeroes in the upper-half complex plane. Since in 
the expression ( lAH both / and g are allowed to have any number of zeroes, this shows that the condition i23\ need not always 
define a good global gauge. 




FIG. 5: Parametric plot of JA^ for «; = 0.1, 0.5, 0.9, 0.999. Increasing K increases the amplitude of h{y), and leads to a singularity at k = 1. 
The function Z{w) has been shifted s.t. h{0) = 0. 
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To describe the deformed surfaces of section llVl we write 

• 1 cos6'o — 1 , . 

9 = .90 + g , with go = — . (A3) 

sm t^o 

The unperturbed planar surface corresponds to g = 0. Other choices of g, which are holomorphic in the upper-half plane 
(including the point at infinity) and for which go + g has no zeroes, describe globally-parameterized deformed fluid surfaces. As 
a simple example, let 6'o = 7r/2 and take 

Y{w) = - f-Kc2™ + log (1 + KC^™)] 

Z{w) = [Ke2™ + log (1 + Ke^™)] 

where k is a real parameter between and 1, and in the expressions for Y and Z we have dropped an irrelevant constant [which 
can be absorbed in a redefinition of the origin of coordinates]. For small k, this function describes a periodic minimal surface 
with period Ay ~ 2tt, and with a deformed contact-line given by h{y) — Kcosy + 0{k^). For k finite, the contact-line profile 
is a complicated function given implicitly by eqs. ( IA4t . and plotted on figure|5l Inserting the above Y and Z in the expression 
( l40l l for the energy gives: 

i?/period = ^ [k^ _ log(l - . (A5) 

This reduces to the JdG energy at small k, and can also be verified numerically. Note that when k ^ 1 the surface becomes 
singular, and the energy per period diverges. 



APPENDIX B: MORE GENERAL LARGE- VOLUME CUTOFF 



In this appendix we will repeat the calculation of the energy of section HVl using a more general container with an outer wall 
at an arbitrary inclination angle. The characteristic function 8c (f^ now reads 



ec(r) = e{x)e[y+^ 



f-^e(^^-y^e{L,~xcoa(l) + zsm(j>) . (Bl) 



The inclination angle cj) of the outer wall is a control parameter, which should drop out in the L^jLy — > oo limit. The contact 
angle of the planar surface with this outer wall is equal to vr — — 6*0, so Young's equilibrium condition requires that the 
corresponding solid-fluid tension be 7" — — cos((/) + ^o)- Repeating the same steps as in section irVI leads to the general 
expression for the energy 

E[h] ^ f Af-d^Ar~-f(f> Ar ■ d^fo + j (f> \h ■ d^ro\-^ h ■ Ar + E^^ry , (B2) 

^ JdD JdD JdV 

where T) = [0, tq] x [— Ly/2, Ly/2\ is the parameter domain defined by 8c(n)) = 1, di_ is the derivative in the inward normal 
direction to dV, and n is the three-dimensional vector normal to the container boundary. 

We can now verify that the inclined wall does not contribute to the above expression. This follows from a fine cancellation 
between the three last terms in eq. (IB2> : 



7 cos Wo z + 



7 sm (p i 7" sin (?o 



1(0 + 0o) 



. (B3) 



We here used the normal vector n ~ {— cos </>, 0, sin 0), which implies that \n ■ Sj^ro | = sai{(f) + 6*0), as well as some three- 
dimensional geometry which is required to extract the contribution of fbnry Doing some straightforward trigonometry, and 
using the fact that 7" = cos(0 + 9^)), one can check that the three terms jB3> indeed cancel. This confirms the decoupling of the 
large-volume cutoff, as was announced in section CVl 
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